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We use tunable laser based Angle Resolved Photoemission Spectroscopy to study the electronic 
structure of the multi-band superconductor, MgB 2 . These results form the base line for detailed 
studies of superconductivity in multi-band systems. We find that the magnitude of the supercon¬ 
ducting gap on both a bands follows a BCS-like variation with temperature with Ao ~ 7 meV. 

The value of the gap is isotropic within experimental uncertainty and in agreement with pure a 
s-wave pairing symmetry. We also observe in-gap states confined to ki? of the a band that occur at 
some locations of the sample surface. The energy of this excitation, ~ 3 meV, is inconsistent with 
scattering from the tt band. 
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Multi-band superconductors have attracted renewed 
interest because of the recently discovered iron based high 
temperature superconductors, whose pairing symmetry 
and mechanism are subjects of a lively debate^. Su¬ 
perconductivity in multi-band systems was already con¬ 
sidered more than half a century ago, immediately af¬ 
ter formulation of the BCS theory, in light of super¬ 
conductivity of transition metal^. It became apparent 
early on that in such a system, the magnitude of the 
order parameter may vary for different bands. On the 
experimental side, beside iron based superconductors, 
several compounds are clearly identified as multi-band 
superconductors, including NbSeJ^, YNi 2 B 2 CP, AsCeJ^ 
and, perhaps most clearly and cleanly, MgB 2 P. MgB 2 is 
one of the most studied multi-band superconductor be¬ 
cause of its relatively high transition temperature (T^ 
^ 40K) and promise as an applied superconducto]!^. 
The superconducting properties of MgB 2 agree well with 
BCS/Eliashberg theory^^i^, and it provides an ideal play¬ 
ground for both theoretical and experimental studies. 
Since its discovery, MgB 2 become regarded as “the” pro¬ 
totypical multi-band BCS superconductor. 

The crystal structure of MgB 2 is very simple, with 
alternating graphite-like boron and magnesium layer^. 
Four electronic bands that cross the Fermi level (Ep) 
can be categorized into two groups: two quasi 2D a 
bands around T and two 3D tt bands near Brillouin 
zone boundar}ffIl. From theory point of view, MgB 2 
is usually described as two band (a and tt) super¬ 
conductor. The pairing is caused by coupling of the 
electron to F 2 ^ phonon mode. However the coupling 
strength is quite different for the two types of bands, 
giving rise to different magnitudes of the superconduct¬ 
ing gap^^^J^. The presence of two different values 
of the superconducting gap was confirmed by several 
spectroscopic measurements^^ Previous Angle Re¬ 
solved Photoemission Spectroscopy (ARPES) measure¬ 
ments have revealed that the gap size on a band is ^6 
meV and the one on tt band is ^3 me\fflf^, which is 
roughly consistent with theoretical calculations. However 


these experiments only measured the magnitude of the 
supercondcuting gap at a single K point on each Fermi 
surface (FS) sheet. No information on how the gap varies 
along one FS has been reported. Electronic Raman mea¬ 
surements with two polarization suggests that the gap 
anisotropy is only 0.4 - 0.6 meV along a F^^. Such 
information is important to understand role of impuri¬ 
ties and interband scattering in prototypical multiband 
superconductor. 

In this paper, we use ultra high resolution tunable laser 
based ARPES^ to study superconducting gap proper¬ 
ties of MgB 2 . Due to limits on the accessible parts of 
the Brillouin zone, the tt band cannot be measured using 
low photon energies. We concentrate our measurements 
on two a bands around F point. Both the temperature 
and momentum dependence of the gap structure are sys¬ 
tematically studied. Our results show that the gap size 
on both a FS sheets are nearly isotropic, consistent with 
theoretical predictions and directly illustrate s-wave pair¬ 
ing symmetry in MgB 2 . Surprisingly, we also discovered 
electronic excitation inside superconducting gap below 
Tc at some locations on the sample surface. 

MgB 2 single crystals with Tc = 39 K were grown by 
a high pressure synthesis technique similar to that de¬ 
scribed in Ref. 23. The typical size of the samples used 
in our measurements is ^0.5 x 0.5 x 0.3 mm^. Samples 
were cleaved in situ at base pressure lower than 8 x 10“^^ 
Torr. ARPES measurements were carried out using a 
laboratory-based system consisting of a Scienta R8000 
electron analyzer and tunable VUV laser light sourc^^. 
All data were acquired using photon energy of 6.7 eV, 
corresponding to = IM12'k jc (inner potential vq ^13 
eV is estimated from ref.l^. The energy resolution of 
the analyzer was set at 1 meV and angular resolution 
was 0.13° and ^0.5° along and perpendicular to the di¬ 
rection of the analyzer slits, respectively. Samples were 
cooled using a closed cycle He-refrigerator and the sample 
temperature was measured using a silicon-diode sensor 
mounted on the sample holder. The energy correspond¬ 
ing to the chemical potential was determined from the 
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FIG. 1. (Color online) Temperature dependence of the superconducting gap in MgB 2 . (a) Sketched FS topology from calcula¬ 
tion. Insert shows the measured FS intensity map close to Brillouin zone center, (b)-(d) ARPES intensity divided by the Fermi 
function along cutl (F-M direction) at 18K, 27K and 40K respectively, (e) ARPES intensity divided by the Eermi function 
along cut2 measured at 22K. Locations of cuts are illustrated in (a), (f) Temperature dependence of EDO’s at of cutl. (g) 
Symmetrized EDCs from panel (f). (h)-(i) Symmetrized EDCs ai kp of cut2 for al and cr2 bands respectively (j) Temperature 
dependence of the superconducting gap extracted from (g)-(f) color/shape coding shown in panel (a). Solid lines indicate BCS 
prediction with corresponding Aq. 


Fermi edge of a polycrystalline An reference in electrical 
contact with the sample. The absence of aging effects 
was verified by thermal cycling. The consistency of the 
data was confirmed by measuring several samples. 

The diagram of the Fermi surface and intensity plot 
at Ej are shown in Fig. la. The ARPES intensity di¬ 
vided by Fermi function plots for several temperatures 
and two cuts are shown in Figs. Ib-e. In the data along 
cutl (Figs. Ib-d), there is only one hole-like band visible 
that corresponds to a2 crossing. Intensity of the inner 
al band is almost completely suppressed due to matrix 
elements. Along cut2 however (Fig. le), both a bands 
are visible, and the the intensity of the inner al band 
is much stronger than the a2. From the location of the 
Fermi crossing point (k^?), we calculate that the radii 
of two a FS sheets are ^ 0.2 and ^ 0.25 re¬ 
spectively. If we ignore the small warping of these two 
sheets along Kz, we can estimate carrier concentration at 
0.069 hole/cell for inner FS sheet and 0.108 hole/cell for 
outer one, consistent with previous quantum oscillation 
result^^^^. 

We next turn to investigating the magnitude of the 
superconducting order parameter on the two a bands. 
The data in Fig. Id was measured above Tc and was 


divided by the Fermi function. It shows a single band 
nicely crossing the Ef characteristic of a metallic state. 
The data in Figs, lb, c and e were acquired below Tc 
and shows the absence of intensity close to Ej and in¬ 
creased intensity at energy of ^ ±7 meV which is the 
hallmark signature of the opening of a superconducting 
gap. The intensity at 7 meV above the Ej is enhanced 
due to thermal excitations above 2A, which is clearly 
visible especially at higher temperatures (Fig. Ic). The 
EDCs at the outer a band are shown in Eig. If. When 
the sample is cooled below Tc, a sharp coherent peak 
emerges at the energy equal to the value of the super¬ 
conducting gap and intensity close to Ej is suppressed. 
Above Ei? a smaller peak is observed due to thermal exci¬ 
tations above 2A. These features are consistent with for¬ 
mation of Bogoliubov quasiparticle^^. The observation 
of the upper branch of Bogoliubov quasiparticle band is 
very useful to double check the location of the chemi¬ 
cal potential which should be centered between the two 
peaks as shown in Eig. If. In order to quantify the 
value of the superconducting gap we symmetrize EDCs 
at kj? (Eig. Ig-i) and fit them with the minimal BCS 
modeP^. Using this procedure we determine the gap 
magnitude at three kp points on the two ES sheets. Po- 
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sitions of each points in Brillouin Zone are marked in 
Fig. la. The resulting gap sizes for each temperature are 
shown in Fig. Ij. All three superconducting gap mag¬ 
nitudes follow BCS-like temperature dependences very 
well with Aq ~ 7±0.5 meV, w hich is slightly larger than 
previous ARPES result^^^^, but more consistent with 
STM resultJ^ and theory calculatiorP^. The gap-to- 
temperature ratio 2Ao/K^Tc' is 4.2, slightly larger than 
the standard single-band weak coupling BCS value 3.5, 
which is likely due to presence of multiple bands. 




FIG. 2. (Color online) Momentum dependence of the super¬ 
conducting gap in MgB 2 . (a)-(b) EDCs at k^ along two a 
FS sheets. The 19 momentum locations in Brillouin Zone are 
marked in left-bottom insert. Data were measured at 25K. 
(c) Angle dependence of the superconduting gap. (d) Same 
data as (c) plotted in polar coordinates after being six-fold 
symmetrized. 

Previous ARPES measurements revealed that the a 
and TT ES sheets have different magnitudes of the order 
parametei!^. But it is not known if and how this quantity 
varies around the respective Eermi surface sheets. To in¬ 
vestigate this we measure the EDCs at k^? around each of 
the (7 sheets and plot them in Eig. 2a, b. We also extract 
the value of the SC gap following the same procedure as 
explained above. Measurements were performed at 25K 
to take advantage of both upper and lower Bogoliubov 


quasiparticle branches. The resulting values of the super¬ 
conducting gap are shown in Eig. 2c, d. The gap mag¬ 
nitudes along both a ES are nearly isotropic and equal 
to ^6 meV (for T=25K). We did not observe any sys¬ 
tematic changes in this value within our ±0.5 meV error 
bars. This result provides direct evidence that MgB 2 is 
a pure s-wave superconductor. Detailed gap calculation 
shows the average gap size on the inner pocket should be 
about 0.5 meV larger than that on the outer pockelP^. 
However at present our experimental accuracy is not suf¬ 
ficient to verify this prediction. Also, since both a bands 
have the same orbital character, the gap difference may 
be possibly suppressed by scattering between these two 
ES sheetP. 



FIG. 3. (Golor online) Inside gap state in superconducting 
state, (a) Measured band structure of cutl( as illustrated 
in fig. la) at 18K. (b) Expanded area close to Ef- Golor 
scale is adjusted to highlight structure above Ef- (c) EDGs 
corresponding to data in (b). Bogoliubov quasiparticle peaks 
and inside gap structures are marked with red and blue arrows 
respectively. 


In fig. 3a, we plot the measured electronic structure 
of cutl in fig.l without removing Eermi-Dirac function. 
A kink structure is clearly revealed at ^ 70meV which 
is caused by electron coupling to E 2 ^ phonon mode. Its 
properties were discussed in great detail previousl}®]. In 
order to reveal the fine structure near Ef, we expand the 
image and adjust the color scale in fig. 3b. Surprisingly, 
in addition to two branches of Bogoliubov quasiparticle 
at ± 7 meV (marked with red arrows), two additional 
intensity peaks exists at ± 3 meV (marked with blue 
arrows). These are states that exist within the larger 
superconducting gap. (Here, the peak positions are de¬ 
termined by local maximum intensity. The same method 
is also used in Eig. 4b.) In order to get more information 
about this in-gap structure, we show the temperature de¬ 
pendent data in fig.4. The color scale in fig. 4 is expanded 
to reveal this very weak structure. As temperature is in¬ 
creased, there are several noticeable changes. The most 
obvious one is that the position of two Bogoliubov quasi¬ 
particle branches moves closer to Ef and the intensity of 
the Bogoliubov quasiparticle branch for positive binding 
energy is increasing. This is because the gap is closing 
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FIG. 4. (Color online) Temperature dependence of the in-gap state. (al)-(a5) ARPES intensity measured at cut 1 for various 
temperatures, (b) EDCs slightly off kp (location marked with red line in al) for various temperatures along line marked in (al). 
Inset shows htting of the EDO above Ef with two gaussian peaks to extract value of the two gaps, (c) Extracted temperature 
dependence of the peak positions of Bogoliubov quasiparticle (red dots) and inside gap structures (blue dots). Two dashed 
lines are guide to the eye. 


and more thermal excitation occur at high temperature. 
At the same time, the position of the inside gap structure 
is also approaching E^. But its intensity becomes weaker 
with increasing temperature and becomes almost invisi¬ 
ble around T^. In fig. 4b, we plot the EDCs slightly off 
kp (location marked with red line in Fig. 4al) at different 
temperatures to quantify this interesting behavior. We 
move slightly off k^ because the strong intensity of the 
main Bogoliubov quasiparticle peak would bury the weak 
inside gap state. Extracted temperature dependent posi¬ 
tions of the Bogoliubov quasiparticle peak and inside gap 
structure are show in fig. 4c. Both of the peak positions 
roughly follow a BCS-like temperature dependence. The 
disappearance of this inside gap structure indicates it is 
closely related to superconductivity. Before discussing 
the possible origin of this new structure, let us summa¬ 
rize its main properties. This structure consists of two 
flat bands at cjq ^ i 3 meV at low temperature, giving 
the ratio of cjo/Aq ^ 0.43. Its intensity is mostly con¬ 
fined at kp in momentum space and its binding energy 
roughly follows BCS-like temperature dependence. We 


emphasize that this in-gap structure is only observed at 
some locations at the sample surface. So far we have 
not found any correlation between this in gap state and 
optically observed imperfections of the surface. 

The presence of this structure is unlikely a trivial ef¬ 
fect of momentum mixing (e. g. due to scattering of 
the photoelectrons) because it occurs only in the prox¬ 
imity of kp. There are several possible explanations for 
presence of the inside gap structure in MgB 2 . The most 
obvious one is inter-band scattering from the tt band due 
to disorder or impurities. The only problem is that pre¬ 
viously measured value of the superconducting gap at tt 
band is smaller by a factor of two i. e. ^1.5 me\®. 
A second possibility is that this structure is an impurity 
bound state in the superconducting gap. Early theoreti¬ 
cal studies of the impurity effect on superconductivity in¬ 
dicate that, at proper scattering conditions between con¬ 
duction electron and impurity, some bound state will ap¬ 
pear inside superconducting gapP®^. The properties of 
this bound state have been intensively investigated both 
theoretically and experimentally because they could give 
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some essential information of superconductivity, such as 
pairing symmetry. Bound states near impurities have 
been observed b oth in conventionaP^and unconventional 
super conduct or in STM measurements. For a con¬ 

ventional superconductor like MgB 2 discussed here, a 
nonmagnetic impurity could not affect superconductiv¬ 
ity much as predicted by early Anderson’s theorenP^, 
which means no bound state would be introduced by a 
nonmagnetic impurity. On the other hand, magnetic im¬ 
purity can give rise to pairs of bound states b oth above 
and below E^?, also called Shiba stateP^^^^^^. The en¬ 
ergy positions of the Shiba states vary with scattering 
details. Up to now, no such states have been reported in 
doped MgB 2 . The Shiba state would be buried under the 
Bogoliubov quasiparticle peak in spectroscopy measure¬ 
ments without momentum resolution. Electronic prop¬ 
erties of the impurity bound state are studied in mo¬ 
mentum integrated density of state in most of the ex¬ 
isting literatures. Information on intensity distribution 
in momentum space is scarce. Whereas this scenario 
seems consistent with our results, one of the key unsolved 
problems is what plays the role of magnetic impurity. 
Both magnesium and boron atoms are nonmagnetic and 
no other magnetic element is involved in sample growth 


procedur^^. Eurther experimental results and theoreti¬ 
cal calculations are needed to identify the scattering po¬ 
tential that plays the role of magnetic impurity. 

In summary, the temperature and momentum depen¬ 
dent superconducting gaps, on the two a bands, were 
systematically studied. The gap size follows a BCS- 
like temperature dependence with Aq ^ 7meV. The mo¬ 
mentum dependent gap structure, on both ES sheets, is 
isotropic, giving direct evidence of s-wave pairing sym¬ 
metry in MgB 2 . We observed a flat-band-hke inside gap 
structure below Tc- Its intensity is mostly confined close 
to kp in momentum space. The energy position of this 
structure also roughly follows a BCS-like temperature de¬ 
pendence. We proposed two possibilities to explain this 
new in-gap structure. Eurther experimental and theoreti¬ 
cal investigations are needed arrive at a definitive answer. 
The electronic properties of MgB 2 revealed here form a 
basis for understanding of multi band superconductivity 
in conventional and unconventional superconductors. 

Research was supported by the US Department of En¬ 
ergy, Office of Basic Energy Sciences, Division of Materi¬ 
als Sciences and Engineering. Ames Laboratory is oper¬ 
ated for the US Department of Energy by the Iowa State 
University under Contract No. DE-AC02-07CH11358. 

* Corresponding author kaminski@ameslab.gov 


^ D. C. Johnston, Advances in Physics 59, 803 (2010) 

^ G. R. Stewart, Rev. Mod. Phys. 83, 1589 (2011) 

^ H. Suhl, B. T. Matthias, and L. R. Walker, Phys. Rev. 
Lett. 3, 552 (1959) 

Y. Kohsaka, C. Taylor, K. Fujita, A. Schmidt, C. Lupien, 
T. Hanaguri, M. Azuma, M. Takano, H. Eisaki, H. Takagi, 
S. Uchida, and J. C. Davis, Science 315, 1380 (2007). 

^ R. J. Cava, H. W. Zandbergen, B. Batlogg, H. Eisaki, 
H. Takagi, J. J. Krajewski, W. F. Peck, E. M. Gyorgy, 
and S. Uchida, Nature 372, 245 (1994) 

® A. F. Hebard, M. J. Rosseinsky, R. C. Haddon, D. W. 
Murphy, S. H. Glarum, T. T. M. Palstra, A. P. Ramirez, 
and A. R. Kortan, Nature 350, 600 (1991) 

^ J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, 
and J. Akimitsu, Nature 410, 63 (2001) 

® G. W. Canfield, Paul C.and Crabtree, Physics Today 56, 
34 (2003) 

^ S. L. Bud’ko, G. Lapertot, C. Petrovic, C. E. Cunningham, 
N. Anderson, and P. C. Canfield, Phys. Rev. Lett. 86, 1877 
( 2001 ) 

D. G. Hinks, H. Claus, and J. D. Jorgensen, Nature 411, 
457 (2001) 

J. Kortus, 1.1. Mazin, K. D. Belashchenko, V. P. Antropov, 
and L. L. Boyer, Phys. Rev. Lett. 86, 4656 (2001) 

J. M. An and W. E. Pickett, Phys. Rev. Lett. 86, 4366 

( 2001 ) 

A. Y. Liu, 1. 1. Mazin, and J. Kortus, Phys. Rev. Lett. 87, 
087005 (2001) 

H. J. Choi, D. Roundy, H. Sun, M. L. Cohen, and S. G. 
Louie, Nature 418, 758 (2002) 

F. Giubileo, D. Roditchev, W. Sacks, R. Lamy, D. X. 
Thanh, J. Klein, S. Miraglia, D. Fruchart, J. Marcus, and 


P. Monod, Phys. Rev. Lett. 87, 177008 (2001) 

X. K. Chen, M. J. Konst ant inovic, J. C. Irwin, D. D. 
Lawrie, and J. P. Franck, Phys. Rev. Lett. 87, 157002 
( 2001 ) 

S. Tsuda, T. Yokoya, T. Kiss, Y. Takano, K. Togano, 
H. Kito, H. Ihara, and S. Shin, Phys. Rev. Lett. 87, 177006 
( 2001 ) 

H. Uchiyama, K. M. Shen, S. Lee, A. Damascelli, D. H. 
Lu, D. L. Feng, Z.-X. Shen, and S. Tajima, Phys. Rev. 
Lett. 88, 157002 (2002) 

S. Souma, Y. Machida, T. Sato, T. Takahashi, H. Mat- 
sui, S.-C. Wang, H. Ding, A. Kaminski, J. C. Campuzano, 
S. Sasaki, and K. Kadowaki, Nature 423, 65 (2003) 

S. Tsuda, T. Yokoya, Y. Takano, H. Kito, A. Matsushita, 
F. Yin, J. Itoh, H. Harima, and S. Shin, Phys. Rev. Lett. 
91, 127001 (2003) 

J. W. Quilty, S. Lee, A. Yamamoto, and S. Tajima, Phys. 
Rev. Lett. 88, 087001 (2002) 

R. Jiang, D. Mon, Y. Wu, L. Huang, C. D. McMillen, 
J. Kolis, H. G. Giesber, J. J. Egan, and A. Kaminski, 
Review of Scientific Instruments 85, 033902 (2014) 

L. Petaccia, C. Cepek, S. Lizzit, R. Larciprete, R. Macovez, 

M. Sancrotti, and A. Goldoni, New Journal of Physics 8, 
12 (2006) 

E. A. Yelland, J. R. Cooper, A. Carrington, N. E. Hussey, 
P. J. Meeson, S. Lee, A. Yamamoto, and S. Tajima, Phys. 
Rev. Lett. 88, 217002 (2002) 

A. Carrington, P. J. Meeson, J. R. Cooper, L. Bali- 
cas, N. E. Hussey, E. A. Yelland, S. Lee, A. Yamamoto, 

S. Tajima, S. M. Kazakov, and J. Karpinski, Phys. Rev. 
Lett. 91, 037003 (2003) 

H. Matsui, T. Sato, T. Takahashi, S.-C. Wang, H.-B. Yang, 







6 


H. Ding, T. Fujii, T. Watanabe, and A. Matsuda, Phys. 
Rev. Lett. 90, 217002 (2003) 

M. R. Norman, M. Randeria, H. Ding, and J. C. Cam- 
puzano, Phys. Rev. B 57, R11093 (1998) 

M. lavarone, G. Karapetrov, A. E. Koshelev, W. K. Kwok, 

G. W. Grabtree, D. G. Hinks, W. N. Kang, E.-M. Ghoi, 

H. J. Kim, H.-J. Kim, and S. I. Lee, Phys. Rev. Lett. 89, 
187002 (2002) 

I. 1. Mazin, O. K. Andersen, O. Jepsen, O. V. Dol¬ 
gov, J. Kortus, A. A. Golubov, A. B. Kuz’menko, and 
D. van der Marel, Phys. Rev. Lett. 89, 107002 (2002) 

D. Mon, R. Jiang, V. Taufour, R. Elint, S. L. Bud’ko, P. G. 
Ganfield, J. S. Wen, Z. J. Xu, G. Gu, and A. Kaminski, 
Phys. Rev. B 91, 140502 (2015) 

Y. LU, Physics 21, 75 (1965) 

H. Shiba, Progress of Theoretical Physics 40, 435 (1968) 


A. V. Balatsky, 1. Vekhter, and J.-X. Zhu, Rev. Mod. 
Phys. 78, 373 (2006), 

A. Yazdani, B. A. Jones, G. P. Lutz, M. E. Grommie, and 
D. M. Eigler, Science 275, 1767 (1997) 

S. H. Pan, E. W. Hudson, K. M. Lang, H. Eisaki, S. Uchida, 
and J. G. Davis, Nature 403, 746 (2000) 

H. Yang, Z. Wang, D. Eang, Q. Deng, Q.-H. Wang, Y.-Y. 
Xiang, Y. Yang, and H.-H. Wen, Nat Gommun 4, (2013) 
P. W. Anderson, Phys. Rev. Lett. 3, 325 (1959) 

H. Shiba, Progress of Theoretical Physics 50, 50 (1973) 

G. P. Moca, E. Dernier, B. Janko, and G. Zarand, Phys. 
Rev. B 77, 174516 (2008) 

J. Karpinski, N. Zhigadlo, S. Katrych, R. Puzniak, K. Ro- 
gacki, and R. Gonnelli, Physica G: Superconductivity 456, 
3 (2007) 


